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I. 



INTRODUCTION 



Compensation or linear control systems by the classical 
methods cf root locus and frequency response is adequate 
when there is a single variable. However when there are two 
or more variables compensation of physical systems for 
desired system dynamics becomes a complex task [1]. 

In modern control engineering systems characteristic 
equations with two or more unknown parameters result. In 
1959, Hitrovic's Method [2] was developed to portray the 
relationship between the roots of the characteric equation 
and two coefficients of the system. This method allowed the 
two lowest ordered coefficients of the characteristic 

equation to be varied. Siljak expanded this method to tae 
Coefficient Plane to allow the variation of any two of the 
coefficients of the characteristic equation. For example: 

n 

i j k 

f (s) = Bs + 3s + ) as (1-1) 

i j l— i k 

k=0 
k*i, j 

where: a are real constants 

k 

B and B are real variables 
i 3 

Application of this method resulted in constant zeta, omega, 
and sigma curves on the B -B plane specifying the roots of 

i j 

equation (1-1) for any choice of £ and oj . 

n 

In 1964, Siljak [3] extended Mitrovic's work so that 
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the two variable paraa.eters may appear linearly in the 
coefficients of the characteristic equation, consider the 
linear case of the characteristic equation: 

n 

f( S ) = ^ a k S = 0 ( T ' 2 > 

k=0 

where: a = bGL + c G + d (1-3) 

k k k k 

and: and are real variable parameters 

Application of this method resulted in constant zeta, omega 
and sigma curves for any choice of Q and ^ on the CL- ^ plane 
specifying the roots of equation (1-2) . 

Siljak's Parameter Plane Method provides the designer 
with a simple procedure for factoring a characteristic 
equation and displaying the results in a parameter plane 
diagram. It allows the designer to obtain information about 
system stability and the affect of parameter variations and 
adjustments on stability. In the general case the Parameter 
Plane Method can be applied to any control problem in which 
it is necessary to determine how the variation of one or 
more parameters effect the root location of a system. 

In 1967, Bowie £4] discovered that a limitation existed 

in parameter plane theory. While attempting to solve a 

sixth order characteristic equation, Bowie found that he 

could not find a complete set of roots in an area of the 

parameter plane that should have contained only real roots. 

By substituting the values of CL and ^ into the 

characteristic equation he found the roots. It was 

discovered that a set of complex roots existed corresponding 

to a £ - CO pair in the left-half s-plane. Further 
n 
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investigation showed that this £-60 pair formed a straight 

n 

line in the parameter plane. This was contrary to the then 

existing parameter plane theory which stated that a £-6 0 

n 

pair will only occupy a point in the parameter plane. 

Bowie called this line of constant £ -constants the 

n 

SINGULAR LINE as it was formed by a singular value of 

£ and GO The singular line added a new dimension to the 
n 

parameter plane because it allows the designer to hold one 
root location constant while moving others to achieve 
desired system characteristics by varying & and A, 
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II. THE SINGULAR LINE IN THE PARAMETER PLANE 



A. THE PARAMETER PLANE 



The Parameter Plane Method was developed by Siljak [3] 
as an extension to Mitrovic's Coefficient Plane. The 
parameter plane is a mapping of points in the s-plane into 
the G< ~ (£> plane. The following is the mathematical 
development of that transformation. 

Consider the characteristic equation: 



where : 
and : 




(2-1) 



a (k-0 , 1 , 2 , . . . / n) are real variables 
k 



s 



— <7+ j (0 





1 - 



2 



( 2 - 2 ) 



Examination of figure 2.1 shows 



and 



-£ = Cose 




Sin e 



(2-3) 



By substituting equation (2-3) , equation (2-2) may be 
expressed as: 



s = CO 


(Cos 6 + 


jSinfl) 


(2-4) 


n 








k k 








and s = CO 

D 


(Cosk Q + 


■jSinkfl ) 


(2-5) 
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+ Im 




Figure 2.1 

A POINT IN THE S-PLANE 
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Since 



it follows that 



and 



-1 

<9 = Cos (-£) 

Cosk 0 = Cosk[Cos(-£) ] 
Sink a = Sink[Cos(-£) ] 



( 2 - 6 ) 



(2-7) 

(2-8) 



Equation (2-7) can be expressed as the Chebyshev function of 

the first kind, T (-U and similarly equation (2-8) can be 

k 

expressed in terras of the Chebshev function of the second 

kind. Sin U £- ^ ]. Substituting equations (2-7) and (2-8) 
k 

equation (2-5) becomes 



k k 

s = CO [T (-O + jSin U (-£) ] (2-9) 

n k k 



k 

where: T (-(){= (-1J T U) 

k k 



( 2 - 10 ) 



and 



u (-n = 

k 



Substituting equation 
substituting equation 
(2-1) becomes 



k+1 

(- 1 ) 0 , (£> 

k 

(2-10) into equation (2-9), then 
(2-9) into equation (2-1), equation 



. n 


n 




k+ 1 

1) 0 (J) 

k 




V -1 k k 

f (s)= ) a CO (-1) I (J) + jsin 
Z_j k n k 

k=0 


E* 

k=0 


k 

CO (- 
k n 


= 0 (2-11) 


for any non-zero value of 


both 


the 


real and 


imaginary 


parts of equation (2-11) 
therefore. 


must 


go 


to zero independently 
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0 



£ 



a 

k 



k (-1) k T (I) = 

n k 



Z i k+ 1 

a „(-1) 0 (1) = 0 

i k k k 



(2t12) 



The Chebyshev functions of the first and second kind 
are given by the recursion formulae 

T (X) = 2T (X) - T m 
k+1 k J k-^ 

(2-13) 

U (X) = 2U( ) - U (i) 
k+1 k-1 



where: 


T (*) = 1 


o m 




II 


0(1) 

0 






0 (3) 
1 



To eliminate T (J) from equation (2-12) it is noted that 
k 



V*» = \ m ’ Vi ro 



(2-K) 



then equation (2-12) may be rewritten 



n n 

v~> k Y. C~» k k k-1 

) (-1) aWC U (X) ] - ) (-1) a, WO (J) =0 (2-15a) 

«L_» k n k k n 



k=0 



k=0 



n 



S" 



k+1 k 
1) a WO O) 
k n k 



= 0 



(2- 15b) 



Examination of equation (2-1 5a) shows that its first part is 
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a linearly dependent solution to equation (2-15b) for ail 
non-zero values of zeta therefore equation (2-15) becomes 



t 

k=0 



n 



k k 

M> a v U V 1 = 0 

k n k-1 



E k k 

(-1) a 0 (p = 0 

k n k 



(2-16) 



k=0 



For the linear case: 



a = db + fic + d 
k k r k k 



(2-17) 



where: 
and : 



b , c , and d are real constants 
k k k 

CL-and ^ are the variable parameters 



Substitution of equation (2-17) into equation (2-16) yields: 



B CL- + C G> + D 
1 V 1 



= 0 



(2-13) 



and 



B a + C ft + D 
2 2 V 2 



= 0 



wher e : 



B 






E k k 

(-D b u m 

k n k 



k=0 .. 



k n k- 1 



k=0 



E k k 

(-1) c s U v Jp 



k=0 



k n k-1 



n 

£ k k 

(-1) c U (X) (2-19) 

^ k n k 
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n 



D 

1 



n 



E k k 

(-1) d 0 



k=0 



k n k-'l ^ 



D 

/ 2 




k 



n 




Application of Cramer’s Rule for the solution of 
simultanous linear equations to equation (2-18) yields the 
desired parameter plane solution equation: 



dt = 



CD -CD 

12 2 1 

6 C - 3 C 
12 2 1 




B D 
2 1 



B C 

1 2 



B D 

1 2 



B C 

2 1 



( 2 - 20 ) 



Siljak [3] defined the parameter plane as a rectangular 

plot with as the abscissa and as the ordinate. Equation 

(2-20) gives CL and A as a function of 1 0 and X . By fixino UJ=a* 

> n i 

and varying (J over the range from zero to infinity in 
n 

equation (2-20), a curve .is produced in the c<-£ plane for 

the fixed value of zeta. This curve specifies the and£ 

pairs which will cause equation (2-1) to have a pair of 

roots corresponding to the fixed value of zeta. Similarly 

by fixing CO ~ CO and varying zeta from -1 to +1 in 

n ni 

equation (2-20), a curve is produced in the CL - ^ plane for 

the fixed value of cj • This curve specifies the CL and $ 

ni 

pairs which cause equation (2-1) to have a pair of roots 

corresponding to the fixed value of ^ . Thus equation 

n 

(2-20) allows the mapping of complex points in the s-plane 
to the parameter plane. 

To transform the real axis points in the s-plane into 
the parameter plane, replaces s in equation (2-1) by to 
obtain : 
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k 



0 



( 2 - 21 ) 



f (s) 



n 



E 



a (-a) 

k 



Now equation (2-17) is substituted into equation (2-21) to 
obtain: 



E 

Ir— n 



( bd+pc + 

k r k 



k 

d) (-a) 



0 



Simplifying equation (2-22) ootains: 



( 2 - 22 ) 



where: 



B(<t) +pC(<r) + D(<r) = 0 



n 



B(a) 



■ E 



k=0 

n 



C(a) 



■ E 



k=0 



k k 
(-1) b a 

A 



k k 

(-1) C a 

k 



D(a) 
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■ E 



k=0 



k k 
(-D a * 
k 



(2-23) 



(2-24) 



For a given value of a the functions 3(a) , C(a) , and D(a) 
are constants and equation (2-20) results in a straight line 
in the parameter plane with a locus of point corresponding 
to the real roots s=-a. 



A digital computer program PARAM A was written by 
Nutting [5] to aid in solving the characteristic equation 
and plotting the constant zeta, omega, found at the end of 
this thesis and a brief explanation is found in Appendix A. 



19 



EXAMPLE I: 



Consider a system 


whose characteristic equation is the 


fourth order polynomial 


• 








4 3 2 


+ (12.5(1+25^ + 42. 








f(s)=s+8.5s +(5CL+28)s 


5) s + 50^ + 25 


=0 


2-25) 


Figure 2.2 shows the 


parameter plane for the 


following 


values of zeta, omega. 


ana sigma: 








OMEGA 


ZETA 


SIGMA 






3.0 


0.00 


0.5 






4.0 


0.30 


1.0 






5.0 


0.35 . 


2.0 






6.0 


0.40 


3.0 






7.0 


0.45 


4.0 






8.0 


0.50 


5.0 






9.0 


1.00 


6.0 






Examination of figure 


2.2 shows the region 


of 


stability 


bounded by ^=0.0, '$ = 1. 


0, and (T“0.0. 


All values 


of dC and ^ 


inside this region 


will result in 


a stable 


system. The 


roots of the system can 


be determined 


by selecting 


a point 


in the plane and reading the corresponding values 


of zeta. 


omega, and sigma. Zeta 


and omega will 


form a 


complex of 


roots so that: 












/ 


2 






s 


= “Jtd ± jCO V 


)-5 








n n v 








and the sigma point will form the real 


roots of 


the 


form: 



S = ~ o- 



Althougli it cay not be immediatly obvious, there are 
two sigma points through each point in the region. Careful 
interpolation or the plotting of addit ionalines is required 
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to obtain the exact roots of the system. But this is 
generally unnecessary as the designer will interested in the 
general location of a specific pair of roots. 

Four points have been selected from figure 2.2. the 
position of these points was measured with an engineer's 
scale and the resulting value of and were substituted in 
equation (2-25) . The roots of the equation were then solved 
by the use of a polynomial rootfinder program on tne IBM 
360/67. Table II- 1 shows the coordinates of the points, the 
resulting roots and the corresponding values cf zeta, omega, 
and sigma. As can be seen the parameter plane provides a 
simple means of determining a root location while variable 
parameters are adjusted. 



TABLE I 1-1 



POINT 


ALPHA 


BETA 


COMPLEX ROOTS 


REAL ROOTS 


OMEGA 


ZETA 


A 


6. 1 


4.0 


-1.461±j5.814 


-1.557 

-4.020 


5. 995 


.244 


B 


6. 1 


1 

3. 0 


-1.615±j5.783 


-1.076 
-3. 794 


6.061 


.299 


C 


7. 0 


1.975 


-2. 278± j6. 075 


-0. 998 
-2.945 


6.490 


.351 


D 


2. 8 


0.725 


-2. 308± j3. 982 


-1.004 

-2.861 


4.603 


.501 



POINTS FROM FIGURE 2.2 
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-scale: 4 units/inch 
-scale: 1 unit/inch 

Figure 2.2 

PARAMETER PLANE EXAMPLE 
FOURTH ORDER POLYNOMIAL 
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THE SINGULAR LINE 



As stated in Chapter I Bowie [4] discovered tne 
existance of the singular line in the parameter plane. In 
investagating the cause of singular lines Bowie found that 
all previous work in the parameter plane was correct for as 
far as it went. It did not, however, give a full 
explanation of the case when both the numerator and 
denominator of equation (2-20) formed singular matrices. 
Siljak [3] did address the fact that this condition could 
exist, nut the idea of a singular line in the parameter 
plane had not, up to that time, been conceived. 

Bowie found that in order for singular lines to exist 
in the parameter plane the following conditions must be met: 



where: 



3 C -PC =0 

12 2 1 



CD -CD =0 
12 2 1 



B D - B D =0 
2 1 12 



(2-26) 



n 



E k k r — i k k 

(-1) b . o nr ) B =) (-1) b, u (Y ) 

k ns k- 1 J s 2 l—i k ns k J s 



B = 



k=0 

n 



E jC iC r — ^ 

(- 1 ) V 0 (X ) C =) (-1) c u (f ) (2-27) 

k ns k- 1 J s 2 i — * k ns k * s 



k=0 

n 



k=0 

n 



Z k k r— i k k 

(-1) d U 0 > D 0 =) (- 1 ) d « ) 

k ns k-1 s 2 Z._/ k ns k s 



k=0 



k =0 
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ind: CO and X are the constantw-constantY values, 

ns s 



equations (2-18) are solved for in terms of £ the 
It is two linearly independent equations: 




Equation (2-28) is the singular line solution. It 
s a straight line in the parameter plane with a slope of 

C and ix intercepts the -axis at -D /C . All points on 

1 1 1 

line represent a singular value of zeta and a singular 
e of omega and therefore map into a single pair of 
lex points in the s-plane where: 



s = 00 ± 

s ns 




(2-29) 



It is obvious that singular lines do not exist for all 

acteristic equations and when they do exist they will 

;t only for the -0U pairs which will satisfy equation 

s ns 

’6) . To aid in determining the exisxance of singular 
;s and plotting them in the parameter plane the computer 
jram SINGULAR LINE has been written. 

Briefly, the program uses the real constants of the 

fficients of the characteristic equation and the desired 

ue of V to solve for value for any real positive values 
3 s 

0 in equation (2-26). The Y and CO values are then 
ns s ns 

<d to generate equations (2-27) and these results are used 

test for the existance of singular lines. If singular 
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region there will be two real roots and a pair of complex 
roots corresponding to the constant ze ta-constant omega 
values of the singular lines. 

A pair of points aas been selected along two of the 
singular lines. Points A and B are on the line for 5 =0.7, 

60=13.28 and point C and D are on the line for J = 0.4 
n 

60=12.309. Table II-2 shows the coordinates, the resulting 
n 

roots and the values of zeta and omega. For Table II-2 the 
coordinates of the point were measured with an engineer’s 
scale and the values of alpha and beta were substituted into 
equation (2-30) and the roots were computed using a 
polynomial rcotfinder subroutine on the IBM 360/67. 

Table II-3 is for the same points, the only difference 
being that the coordinates were computed using the 
slope-intercept method utilizing the printed output from the 
program SINGULAR LINE. 

Table II-3 shows that this method of obtaining system 
roots provides an accurate method for obtaining singular 
lines and that the lines a for constant zeta-consta at omega 
values, while Table I I- 2 shows that reasonably accurate 
results can be obtained using a minimum of equipment. 
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“-scale: 5 units/inch 
/?-scale: 1 unit/inch 

Figure 2.3 

SINGULAR LINES CN THE PARAMETER PLANE 
FOURTH ORDER POLYNOMIAL 
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TABLE I 1-2 



POINT 


ALPHA 


BETA 


COMPLEX ROOTS 


REAL ROOTS 


OMEGA 


ZETA 


A 


20. 


19.67 


-9.244± j 10. 87 


-.9326 
- 422.6 


13.40 


. 639 


B 


10. 


10.5 


-9.281±j9.670 


-1.051 

-222.4 


13.40 


.693 


C 


10. 


5.667 


-4 . 83 9± j 1 1 . 1 9 


-.9971 
-231 . 3 


12.20 


. 39 7 


D 


IS. 


8. 167 


-4 . 9 1S± j 1 1 . 27 


-1. 045 

-326.7 


12.30 


.400 



POINTS FROM FIGURE 2.3 



TABLE II-3 



POINT 


ALPHA 


BETA 


COMPLEX ROOTS 


REAL ROOTS 


OMEGA 


ZETA 


A 


20. 


13-66 


-9 . 298± j9. 486 


-1.071 

-422.3 


13.40 


.700 


B 


10. 


10.52 


-9.298±j9. 486 


-1.051 

-222.4 


13.40 


.700 


C 


10. 


5.755 


-4. 920± j 11.27 


-.9982 

- 23 1 . 2 


12. 30 


.4 00 


D 


15. 


8. 225 


-4 . 920± j 1 1 . 27 


-1.011 

-331.2 



12.30 


.4 00 



POINTS FROM FIGURE 2.3 
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Ill 



• PR OPE RTIES OF THE SINGULAR LINE 



To fully understand the relationship of singul 
to the parameter plane it is necessary to examine 
singular line is and what ix represents. 

The singular line exists as a unique situa 
occurs while mapping the s-plane into the paramete 
The necessary condition that must exist is 
determinants of the matrices which form equation (2 
be singular. tfhen this occurs, a discontinuity a 
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not always occur. 



EXAMPLE III 

Consider the following aircraft pitch stabilization 
system [ 6 ]. 




Figure 3.1 

AIRCRAFT PITCH STABILIZATION SYSTEM 

where: $ is the pitch rate command 

i 

0 is the pitch rate 

0 is the pitch angle 
o 

^ is the elevator angle 

and: G ^ (s) is a pole-zero compensator 

G (s) is a hydraulic servo 
2 

G^(3) is the airframe 
H (s) is a rate gyro 
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The transfer functions of the system components are: 
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The characteristic equation for the uncompensated system is: 

5 4 3 2 

s + 62.38s + 374 1.01s + (31092.5 + 61732. 6K)s 





+ (398891 + 1 9877 9K) s = 0 (3-2) 


where K = K K 

1 2 




This system 


is unstable for K>3.27. It is desired to 


compensate for K = 


5 and to avoid the natural frequencies of 


the airframe, qj = 

n 


3.24, and the hydraulic servo, £0 = 7.049 

n 


The characteristic 


equation of the compensated system is: 


6 5 


- 4 



s + (62. 38 +£)s + (3741.01+ + 62.38@)s 

3 

+ (31092.5 + 61732. 8K + 374 l.01^)s + (398891 + 198751K<£ 

2 

+ 31092. 5^) s + (198751 K2C+ 398891£)s = 0 (3-3) 

This characteristic equation with K = 5 was used with 
the program PARAI1 A to obtain the parameter plane curves, 
figure 3.2, and with the program SINGULAR LINE to determine 
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the existauce of singular line. Two singular lines were 
found for each zeta value of 0.1, 0.2, 0.3, 0.4, and 0.5, 
and one singular line was found for each zeta value of 0.6, 
0.7, 0.8, and 0. 9. 

Zeta = 0.6 was chosen for compensation as it has an 
omega of 12.74 radians and would result in a fairly fast 
response without causing resonance at the natural 
frequencies of the airframe and the servo. The singular 
line for zeta = 0.6 is the broken line on the parameter 
plane of figure 3,2. 

Three points were selected on the singular line. The 
points were; A{14. 56,50), 3 (21.13, 150) 8 and 0(27.69 ,250). 
These values of alpha and beta - where put into the 
characteristic -equation and the roots were solved for on the 
IBM 360 with a polynomial rootfinder routine. 

There were two real roots and two pairs of complex 
roots tor each of the points. ^ Contrary to the supposition 
that al.l points along the singular line will produce a pair 
of complex roots corresponding to the value of the 
zeta-omega of that line, in this case such roots were not 
found. The singular line with zeta =0.6 and omega = 12.74 
should have according to prevoius concepts result in a 
complex rcot pair located at 

s = -7. 645 ± jlO. 194 (3-4) 

This was not found to be the case. The complex roots 
that were found did not agree with the singular line roots, 
but they did agree with the parameter plane roots. These 
results are as follows: 
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“-scale: 20 units/inch 
/?-scale: 50 units/inch 



Figure 3.2 

AIRCRAFT PITCH STABILIZATION SYSTEM PARAMETER PLANE 
WITH ZETA=0.6-OMEGA = 12. 74 SINGULAR LINE 
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point 
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ZETA 


OMEGA 


A 


-8.083± j3. 71 8 
-2.512 + 1)69.00 


. 908 
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8.897 

69.04 


B 


-7 . 944± j 6 . 987 
-16. 6 6±j 64 . 4 1 


.751 
. 251 


10.58 

66.53 


C 


-7.871±j8.0l2 
- 20. 75±j 59. 36 


.701 

.330 


11.23 

62. 87 



Identical results were found when values were taken 
from other singular lines in the same relative area of the 
parameter plane. All the roots were found to agree with the 
parameter plane, but not with the singular lines. 

The initial response to these results was that the 
singular lines that were found were incorrect. However, 
after carefully checking the program siat were found yere 
incorrect. However, after carefully checking the program 
SINGULAR LINE and hand calculating two of the singular 
lines, the program was deemed to be correct. 

In an attempt to find a solution to this dilemma the 
singular line examples given by Bowie were examined. In all 
cases che singular lines were located outside the region of 
the normal complex root parameter plane curves. 

In the belief that the constant zeta curves in figure 
3.2 formed a loop as omega increased, values of alpha and 
beta believed to oe outside the region were tried. The 
straight lines connecting the zeta curves in the figures 
shown in this thesis are the result of the DRAW routine and 
are not part of the parameter plane. They result when 
intermediate points are not in the area of the plot. Again 
three points on the same singular line were selected. These 
points were D (76 . 93 , 1 000) , E (142.57,200 0) and 
F (208 . 22 , 300 0) . These points are not located on figure 3.2, 
but the location was computed using the slope- intercept 
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method. These results are as follows: 
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relatively insensitive to large changes in alpha and beta. 



In light of these results an examination of the 
parameter plane in the region of these points was conducted. 
For reasons of clarity only three values of zeta were 
plotted. These were 0.0, 0.6, and 1.0. The resulting 
parameter plane is seen in figure 3.3. As can be seen the 
constant zeta lines do not loop around as prevoiusly 
thought, but are discontinuous and the singular line 
approaches asymptotically. 

What has happened is that as omega was increased on the 
constant zeta curve and approched the singular point the 
matrix which forms the denominator of eguation (2-20) 
becomes very much smaller than the numerator causing the 
constant zeta line to go to plus infinity at the singular 
point asymptotic to the singular line. Once the value of 
omega is past the singular point it reenters the parameter 
plane at minus infinity, asymptotic to the singular line. 
Since there is only one singular point for this value of 
zeta the constant zeta line will continue through the 
parameter plane approaching infinity once more as omega 
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approaches infinity. 

The reason that this occurs is evident when the mapping 

of the constant zota line form the s-plane to the parameter 

plane is observed. The constant zeta line in the s-plane is 

continuous foria CO =0 to CO = oo, therefore the mapping of this 

n n 

line into the parameter plane must also be continuous. 
However there is a discontinuity at the singular point where 
the constant zeta curve in the paramter plane goes to 
infinity for a finite value of omega. The singular line 
provides a continuous mapping of the constant zeta curve as 
it traverses from plus infinity to minus infinity. This 
line exists at infinity, but not in the parameter plane as 
it is known. 

In an effort to examine the change in alpha and beta as 
omega was increased a printout of PARAM A was made. However 
due to the iteration of omega in PARAM A it was not possible 
to obtain sufficient data around the singular point. In 
PARAM A for an omega of about 12 radians the iterative step 
is about 0.9 radians. Therefore the section of PARAM A 
which computes the constant zeta curves was modified to vary 
omega over 1 radian in steps of 0.001 radians starting 0.5 
radians below the singular point and increasing to 0.5 
radians above the singular point. 

In order to see how small increases in omega affect the 
values of alpha and beta some of the points are shown: 
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Figure 3.3 

PARAMETER PLANE FOR AIRCRAFT PITCH STABILIZATION SYSTEM 

ZETA =0.0 
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As can be seen for large changes in alpha and beta in the 
area of the parameter plane where the singular line is 
asymptotic to the constant seta curves results in small 
changes of omega. Since omega is very insensitive to 
changes or: alpha and beta in this region the concept of the 
singular line is valid from the standpoint of the engineer. 
It is possible to vary a parameter along the singular line 
in this region and expect to maintain a constant value of 
zeta and omega. 

EXAMPLE IV 



In light of the unexpected results obtained in example 
III it was decided to carefully reexamine example II that 
was given in the previous chapter. The parameter plane was 
expanded around the area of the origin where the singular 
lines pass through the complex root portion of the parameter 
plane . 

Figure 3.4 shows this area of the parameter plane. The 
constant zeta curves are for values of 0.0, 0.2, 0.5, 0.6, 
0.8, and 1.0. There are two singular lines for each zeta 
value of 0.2, 0.5, 0.6, and 0.8 which are shown as dashed 
lines. There are two discontinuities in each zeta curve. 
Although it is not immediately obvious due to the limited 
area of the plot, all the constant zeta curves originate in 
the second guadrant, and become infinite at the first 
dicontiuity and then reappear at infinity in the fourth 
guadrant. The curves, for values of zeta up to and 
including 0.5 then follow a path just to the left of the 
origin and enter the first quadrant where they again become 
discontinuous and infinite. The zeta curves for values of 
zeta greater than 0.5 come from infinity in the fourth 
guardant and go to the left of the origin where they turn 
into the third guadrant where they again become 
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discontinuous and infinite. The curves for zeta values of 
0.5 and less reappear at infinity in the third quadrant then 
go to infinity as a continuous function in the first 
quadrant. The curves for zeta values greater than 0.5 
reappear at infinity in the first quadrant, sweep toward the 
origin then rack through the first quadrant to infinity 

The region of stability in the parameter plane is 
bounded by the zeta=0.0 curve, therefore the singular lines 
of interest lie in the first quadrant. As can be seen the 
zeta=0.5 and zeta=0.6 singular linos nearly coincide with 
the corresponding zeta curves in the parameter plane and tne 
zeta=0.2 and zeta=0.8 lines are being asymptotically 
approached by the corresponding parameter plane zeta curves. 

As in example III the modified section of PARAb A was 
used to abtain a printout of alpha and beta for a 1 radian 
change in omega as the constant zeta curve passed through 
the singular line value. In this example the constant zeta 
curves in the first quadrant approached their singular value 
very rapidly. 

For purposes of illustration points have been placed on 
the constant zeta curves in figure 3.4 at the value at which 
omega is within 0.5 radians of the singular point which is 
at infinity. These points are: 
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0.571 


1.60 


C 


0.80 


4.754 


6.62 



39 




-scale: 2 units/inch 
-scale: 2 units/inch 



Figure 3.4 

SINGULAR LINE THROUGH THE PARAMETER PLANE 
FOURTH ORDER POLYNOMIAL 
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As car. be readily seen these points are within 5 
percent cf the singular value of omega and if one continues 
along the zeta curves it will be found, that for uractrical 
design purposes, the changes in omega are inconsequential. 



It has been shown that a true singular line in the 
original concept aces not exist. In examples III and IV it 
is shown that, if a designer carefully examines the singular 
lines and their relationship to the constant zeta curves of 
the parameter plane the singular concept can be used for 
practical design applications. 

MMELj; V 

This example is the original system in which Bowie 
found that singular lines existed in the parameter plane. 
The system which was used has a characteristic equation 
which is the following sixth order polynomial: 

6 5 4 3 2 
s + 80s + (20* + 1600)s + 840 c<s + ( 1 600ci + 400$) s 

+ 160C£s * 1600$ = 0 (3-5) 

The parameter plane with constant zeta curves with the 
values of 0.0, 0.2, 0.4, 0.6, 0.3, and 0.9 is seen in 
figures 3.5 and 3.6. Figure 3.6 is the area around the 
origin in figure 3.5. There were two singular lines tor 
each of the zeta values 0.2, 0.4, 0.6, 0.8, and 0.9 and 
these are shown on figures 3.5 and 3.6 as dashed lines. 

As can be seen the parameter plane has no visible 
discont inuities as omega increases, however due to the fact 
that two singular lines exist for each value of zeta there 
roust be two discontinuities in each zeta curve. Examination 
of the singular lines ana the constant zeta curves of the 
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parameter plane indicates chat the singular lines are 
tangent tc the corresponding zeta curve at one point. When 
the values of zeta and omega for these points were examined 
it was found that the tangent points occur at the singular 
points in the constant zeta curves. 

Unlifce examples III and IV where the mapping from the 
s-plane approached infinity as the singular point was 
approached this example maps as though no discontinuity 
existed. When the discontinuity is reached the constant 
zeta curve immediately goes to plus infinity tangent to the 
curve at the point of discontinuity and returns from minus 
infinity tangent to the constant zeta curve and continues 
mapping to infinity as omega approaches infinity. 

This phemonenon does not reveal itself in a normal 
parameter plane study for two reasons. The first of these 
is the previously mentioned size of iteration in the PARAM A 
program and the second reason is that a normal practice in 
computer programing is to limit the numerator of a term to a 
minimum size to prevent an overflow or underflow in the 
computer when a number beyond the computer’s capablity is 
reached. This practice was followed in the writing of PARAM 
A so if the singular point was reached on an iterative step 
the program would have recognized this as a forbidden number 
and continued to the next value of omega, thereby bypassing 
the singular point and continuing with no indication that 
this condition existed. If Bowie had not oeen interested in 
the area outside the normal complex root area of this 
example he would have had no indication that singular lines 
existed. 
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Figure 3.5 

SINGULAR LINE IN THE PARAMETER PLANE 
EXAMPLE V 
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Figure 3.6 

SINGULA F; LINES IN THE PARAMETER PLANE 
EXAMPLE V - EXPANDED 
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Ia this example the singular lines are true singular 
lines as originally conceived. If a value of alpha and beta 
is selected at any point along these lines a pair of complex 
roots corresponding to the constant zeta-constan t omega of 
that line will result. To show this fact six points have 
been selected in the region of singular line on figure 3.5. 
As can be seen by the following the expected complex roots 
were obtained. 
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0.80 

0.40 


23.55 
47.4 1 


-18.44±j14. 13 
-18.96±g43. 15 


-18. 42 + j 14.15 
-18. 95±g43.44 


B 


0.90 

0.60 


20.88 

31.53 


-18.79±j9.101 

-18.92ij25.22 


-18. 79±19. 121 
-18.9 1±3 25. 23 
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0.80 

0.60 


3.396 

31.53 


-2.717±j2.078 
-18 .92±-j25. 22 


-2.72U 12.081 
-18.9 1+325.22 


D 


0. 90 
0.60 


3.831 

31.53 


-3.450+10.728 

-18.92+g25.22 


-3.465±10.729 
-18. 90+j25. 23 
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0.60 

0.40 


2.537 

47.41 


- 1 . 522±'j2 . 030 
-18.96+143.45 


-1.524±12.032 
-18. 94±g43.43 


F 


0.60 

0.80 


2.537 

23.55 


-1 . 522±j2 .03 0 
-18.44±114. 13 


-1.521± j2.029 
-18. 4 7±1 14.15 



In summary, the singular lines result when a constant 
zeta or constant omega curve is mapped from the s-plane to 
the parameter plane and a discontinuity occurs. At this 
discontinuity the parameter plane mapping goes to infinity, 
either plus or minus, and returns from the other side of 
infinity. The singular line simply provides a mapping 
through this region to maintain a continuous mapping from 
the s-plane to the parameter plane. In cases, such as 
examples III and IV, where the mapping first approaches 
infintiy as the singular point is approached the singular 
line concept is a valid engineering aid over a limited range 
of the singular line. However in cases such as example V 
the singular line may be used over its entire range as an 
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excellent design aid. 



IV 



. THE SINGULAR SURFACE IN THE PARAMETER SPACE 
A. PARAMETER SPACE 

The parameter space is an extension of the parameter 
plane when there are three or more variable parameters [7], 
[8]/ [9]. Once the concept of the parameter plane is 

understood it is relatively easy to extend the concept to 
three or more variable parameters. 

The easiest method to generate a three dimensional 
parameter space is to set the third variable parmeter the 
third variable parameter to a known value and then generate 
the two dimensional parameter plane. The process is then 
repeated for additional values of the third parameter. The 
result is a series of parameter planes which will form a 
parameter space. 

EXAMPLE VI 

Consider the following system: 




Figure 4 . 1 

THIRD ORDER SYSTEM WITH DUAL CASCADE COMPENSATION 
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The characteristic equation for this system is : 

5 4 3 2 

s + (33 + ^)s + (9 2 + 32£) s + (60 + 92 £ + KJ s 

+ (2 OK + <*K + 60 ) s + Kc< = 0 (4-1) 

Figure 4.2 shows the parameter space for the constant 

zeta=0.5 and constant omeya=1.7 and the values of K from 10 
to 15. It was constructed by fixing the value of K and 
computing the parameter plane using PARAM A. The curves were 
then plotted onto a single plot to show the three 

dimensions. 

Both the constant zeta curves and the constant omega 
curves form a surface and the interconnection of these 
surfaces forms a curved line through the parameter space. 
Any point on this line will generate a pair of complex roots 
corresponding to the s-plane point of zeta=0.5 ana 

cmega=1.7. For these values of zeta and omega the complex 
roots are: 

s = -0.850 ± j 1 . 47 0 (4-2) 

By selecting three points on this line this is shewn to 
be true. The points are A(0.16, 2.6, 10.0). B{0.347, 3.0, 

12.0), and C(0.52, 3.63, 15.0). These poinfs were taken 
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B. 



SINGULAR SURFACES 



The parameter space is continuous throughout the range 
of the third variable parameter. If for any value of the 
third parameter a discontinuity appears in the parameter 
plane it will also appear in the parameter space as a 
discontinuity. As stated in chapter III this discontinuity 
will result in a singular line. As the third parameter is 
varied the discontinuity will result in singular lines for 
the same value of zeta and omega. The slope and the axis 
intercept points will change slighty so the result is a 
smooth surface representing a warped hypersurface traversing 
through the parameter space. 

Figure 4.3 is a photograph of a singular surface model. 
The singular surface is from the system in example VI. the 
singular values are zeta=0.5 and omega=12.04 and the third 
variable K is varied from 5.0 to 20.0. As can be seen the 
surface is a smooth warped hyperplane. The model was 
constructed by obtaining the singular lines for the same 
singular value at each value of K, transferring them to 
cardboard, cutting them and accentting them with a 
contrasting color. The photography was done by the Photo 
Division, Education Media Department, Naval Postgraduate 
School . 

Figure 4.4 is a drawing of the same singular surface as 
figure 4.3. The difference between the drawing and the 
photograph is that in the photograph the K axis starts at 
K=5 in the foreground continues back through to K=20, while 
in the drawing K=5 is the most distant singular line and the 
singular lines progress forward until the value of K=20 is 
reached. 
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Figure 4.3 
SINGULAR SURFACo 
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elationship of the singular surface to the 
space is exactly the same as that of the singular 
parameter plane. As with the singular line, the 
urface provides continuity at infinity. The 
st be aware of the type of discontinuity which is 
the surface. If the discontinuities occur like 
example III and IV then the principles -of 
ines may only be applied in limited areas of the 
f the singularities occur like those in example 
ny point on the singular surface will give a pair 
roots corresponding to the constant zeta-constant 
e surface. 



53 



V_. COHPENS ATION T££HNI2jJES USING THE SINGULAR LINE 

The major advantage of using singular lines for 
compensation is that as one moves up and down the singular 
line, witnin previously mentioned limits, one pair of 
complex roots remains in the same locations in the s-plane. 
While this pair of roots is stationary other root locations 
can be varied as the location on the singular line is moved 
through other singular lines, real root lines, and other 
constant 2 eta parameter plane curves. 

There many excellent texts available which go into 
great detail on how root location will affect system 
response, overshoot, damping and so forth [1], [10]. Once 
the basic principles on how root locations affect a system 
are understood it is possible to use singular line concepts 
to obtain proper root locations. 
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EXAMPLE VIII 

Cons id er the pitch stabilization system given in 

example III. 




Figure 5.1 

AIRCRAFT PITCH STABILIZATION SYSTEM 

The characteristic equation for the system is: 

6 5 4 

s + {62.33 +^) s + (3741.01+ + 62.33 p) s 

3 

+ (31092.5 + 61732. 8K + 3741.01£)s + (398891 + 198751 FOC 

2 

+ 31092. 5^) s + (193751 Kd+ 39S891£)s = 0 (5-1) 

It was shown in example III that singular lines exist for 
this system, but they are the result of discontinuities at 
infinity. Therefore in order to effectively use the 
singular line concept it is necessary to examine the 
singular lines to determine the areas in which singular 
results could be expected. 

In the case of this problem it was decided to use the 
singular line with zeta=0.6 and oraega=12.74 as a dominant 
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root pair to give fast system response and avoid the natural 
freguenci.es of the servo and the airframe. Figure 5.2 shows 
the singular lines for this system for values of zeta of 
0.4, 0.5, 0.6, 0.7, 0.8, and 0.9. 

By moving up the zeta=0.6, omega=12.74 singular line 
above the region where beta=3000 the value of zeta and omega 
should remain relatively constant. Two points have been 
selected to attempt, compensation, A (273. 87, 4000) and 
B (306.6 9, 4500) . These values have been put into the 
characteristic equation and the system was simulated using a 
step input of 0.2 radians per second. 

Figure 5.3 is the response for the point A values of 
alpha and beta and figure 5.4 is the response for the point 
B values of alpha and beta. As can be seen the system 
response is essentially the same for both points. The 
response from the point A values has a peak overshoot of 
1.08 and the response from the point 3 values has a peak 
overshoot of 1.1. The steady state output is 1.66 degrees 
in both cases. 

This demonstrates one of the main advantages of the 
singular line concept. As two parameters are varied the 
system response remains insensitive to the change. 

If this parameter variation were to take place inside 
the normal complex root region of the parameter plane the 
shift of root locations would result in a change of output 
response. Figure 5.5 is the parameter plane of the normal 
complex root area of the airdraft pitch stabilization 
system. Between point A and point C the linear change in 
beta is 200 and the change in alpha is 13.13. 

There are two pairs of complex roots and two real roots 
in this section of the parameter plane. One of the real 
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roots is cancelled at the origin of the s-plane by the zero 
produced b y the rate gyro in the feedback loop. The other 
real root location is past s=- 75 in these cases and will 
only have a slight affect on the steady state value of the 
output. The complex root locations and their corresponding 
values of zeta and omega are: 



POINT 


COMPLEX ROOTS 


ZETA 


OMEGA 


A 


-8.033±j3.713 

-2.512±j69.00 


.908 

.036 


8.897 

69.04 


B 


-7 . 944 + j 6 . 987 
- 16. 66± ] 64 . 4 1 


.751 

.251 


10.58 

66.53 


C 


-7.871+18.012 
-20. 75±g 59.36 


.701 

.330 


11.23 

62.87 



At point A the dominant root is located at s = -2.512 ± 
j69.00. At point C the dominant root is located at s = 
-7.871 ± j8.012. It can be seen that at the point A values 
the system response will be almost undamped with a peak 
overshoot of about 1.9 and a natural frequency of about 69 
radians, while at the point C values the system will be well 
damped with a natural frequency of about 11 radians. 
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“-scale: 200 units/inch 
/3-scaie: 1000 units/inch 



Figure 5.2 

SINGULAR LINES FOR 
AIRCRAFT PITCH STABILIZATION SYSTEM 
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Figure 5.4 

STEP RESPONSE FOR POINT B 
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a -scale: 20 units/inch 

/?- scale: 50 units/inch 



Figure 5.5 

AIRCRAFT PITCH STABILIZATION SYSTEM PARAMETER PLANE 

\ 
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VI. 



CONCLUSIONS AND H EC ON HEN DAT 10 NS FOR FURT HER STUDY 



It has been shown that some of the original concepts of 
the singular line were in error due to previously unknown 
properties of the singular line and its relationship to the 
parameter plane. The main reason these properties were not 
found earlier was due tmainly to standard computer 
programming practices and limitations. 



Although it has been found that some limitations exist 
when using singular line concepts it has been shown that the 
basic concepts of singular lines can be used by the design 
engineer to easily obtain desired system response when more 
than one variable parameter is present in the system. 



An example has been presented which shows that the 
singular line concept is an excellent method to obtain 
desired system response. The method will allow the designer 
to select either dominant root locations or, if it is 
desired, nondominant root locations. the root locations 
when selected from the singular line region of the parameter 
plane will remain fixed in their location as long as the 
alpha and beta parameters are restricted to the singular 
line. 



For further study it is highly recommended that the 
case where true singular lines exist such as in example V be 
thoughly studied to obtain a full understanding of why the 
singular point does not occur at infinity in the parameter 
plane. Once this is accomplished devise a method to force 
this type of condition at a desired zeta-omega value. If 
this can be accomplished it would allow the designer to 
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li ue 



select a 
conce pts 



root location in the s-plane and use singular 
for easy design and compensation. 
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APPENDIX A 



PARAM A 



This is a modified version of the PARAM A program 
written by Nutting [5J. Originally written in FORTRAN for 
the CDC 1604 Computer, it has been modified to FORTRAN IV 
for the IBM 360/67 at Naval Postgraduate School. It has 
also been modified to allow the user to select the size of 
the output plot and to plot singular lines in the parameter 
plane. 



However PARAM A will not determine the existance of 

singular lines. It is necessary for the user to know the 

value of £ and u of a singular line before this option can 

n 

be used. These values can best be determined by use of the 
program SINGULAR LINE. The user then can either put these 
values into PARAM A or use the values to draw the singular 
line on the PARAM A output plot. 
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APPENDIX B 
SINGULAR LINE 



The purpose of the program SINGULAR LINE is to 
determine the existance of singular lines in the parameter 
plane. I he program is written in FORTRAN IV for the IBM 
360/67 computer and the CALCOMP plotter at Naval Postgrauate 
School. 

The values read into the program are the coefficients 
of the characteristic equation, the values of zeta to be 
examined, and plotting information. The program contains 
three subroutines which are located in the source library of 
the IBM 360/67. These are POLRT which has the ability to 
solve up to a 36 order polynomial, DRAW which outputs the 
plot to the CALCOMP plotter, and BCNV which converts numeric 
symbols to alpha-numeric for use in DRAW. 

The program takes the values of zeta and first 

generates the Chebyshev function of the second kind. It 

then generates the coefficients for the omega terms in the 

polynomial expansion of B C - 3 C =0. This generation is 

12 2 1 

accomplished by using the coefficients of the alpha and beta 
terms in the characteristic equation, the value of zeta and 
the Chebyshev function. The resulting omega coefficients 
are now input into the subroutine POLRT. 

POLRT solves for the roots of the omega polynomial and 
outputs the results to the program. These roots are now 
examined for positive real values. If positive real roots 
are present there is a possibility of singular lines 
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existing for that value of zeta. If there are not positive 
real roots, there are no singular lines. This information 
is printed out and the next value of zeta is examined. 

For the positive real roots B , B , C , C , D , and D 

12 12 1 2 

of equation (2-27) are formed and used to insure that a 

singular line does exist. If a singular line does exist 

then B , C , and D are used to generate tne points along 

1 1 1 

the line by using equation (2-28) . These points are now 
input into the subroutine DRAW and are plotted and labeled. 

The printed output includes the value of zeta and 
omega, the alpha and beta intercept points, the slope, and 
the label. If there are less than two points generated 
within the area of the plot there will be no plot of that 
line. This information is indicated on the printed output 
and the user can either change the plot scales or use the 
output line information to construct a singular line plot. 

This procedure is repeated for each value of zeta, 
plotting all the lines on a single plot. Multiple runs for 
different characteristic equations are possible by using 
additional data decks. Each equation will have its own 
output plot. 
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